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Full distributions of conductance through quantum dots 
with single-mode leads are reported for both broken and un- 
broken time-reversal symmetry. Distributions are nongaus- 
sian and agree well with random matrix theory calculations 
that account for a finite dephasing time, t v , once broadening 
due to finite temperature T is also included. Full distribu- 
tions of the derivatives of conductance with respect to gate 
voltage P(dg/dVg) are also investigated. 
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The remarkable success of random matrix theory 
(RMT) and other noninteracting theories in describing 
the statistics of quantum transport in mesoscopic elec- 
tronic systems is surprising considering that electron- 
electron interactions not accounted for in the theory are 
sizable compared to other energy scales in these systems. 
In the last few years, an essentially complete statistical 
theory of mesoscopic fluctuation and interference effects 
in disordered or chaotic (irregularly shaped) quantum 
dots and wires has been formulated using these meth- 
ods [|l|,|U|] and has recently been extended beyond the 
"universal" regime to include short-trajectory effects j|, 
clarifying the degree to which these theories are indeed 
universal ||. 

But how seriously should a noninteracting theory be 
taken when describing real metallic or semiconductor 
structures? Apparently, this depends on the quantity 
being measured. For instance, the mean and variance 
of mesoscopic conductance fluctuations in open quantum 
dots and disordered wires appear to be in good agree- 
ment with random matrix theory [EJ and sigma model 
calculations [El once temperature and dephasing effects 
are included, whereas the distribution of energies needed 
to add subsequent electrons to a closed dot does not ap- 
pear distributed according to the famous Wigner-Dyson 
law, a basic result of RMT 0. 

In this Letter we carry out a stringent test of statis- 
tical theories of mesoscopic conductance fluctuations by 
measuring the full distribution P(g) of conductance, g, 
through chaotically shaped ballistic quantum dots. Con- 



ductance distributions of quantum dots have previously 
been calculated within RMT JIH, including the effects 
of dephasing |l|, pXi| , pT| , ^2| . The distributions are univer- 
sal for any fully chaotic or disordered dot, sensitive only 
to whether time reversal symmetry is obeyed (/3 = 1) 
or broken (/3 = 2), controlled by adding a magnetic 
field, B > ~ 4>o/Adot, where A^ot is the dot area and 
4> = h/e is the flux quantum. RMT yields quite inter- 
esting (i.e. strongly nongaussian) distributions when one 
or two quantum modes connect the dot to bulk reser- 
voirs. To date, however, experimental measurements of 
these nongaussian distributions have not been reported, 
first because it is difficult to generate large ensembles 
of statistically identical devices, and second because de- 
phasing, which acts roughly as extra modes coupling the 
dot to the environment, leads to nearly gaussian dis- 
tributions nng. To see the nongaussian distribu- 
tions, dephasing rates and temperatures comparable to 
the quantum level spacing are required. We solve the first 
problem, of obtaining large ensembles, by using electro- 
static shape distortion of gate-defined quantum dots in a 
GaAs/AlGaAs heterostructure p3| . 

We find good agreement between the experimental dis- 
tributions and the RMT predictions over a broad range 
of temperatures once thermal averaging is properly ac- 
counted for. The agreement is particularly surprising 
since we are investigating the case of single-mode leads, 
N = 1. This is the transition between open and closed 
dots; for any lower conductance to the dot, electron- 
electron interactions in the form of Coulomb blockade 
dominate transport, leading to dramatic departures from 
a noninterating picture pSj . Here N denotes the num- 
ber of modes or channels in both the left and right leads, 
giving dimcnsionlcss lead conductances g\ = g r — 2N , 
where g is in units of e 2 /h. 

Distributions of mesoscopic conductance at T — may 
be calculated within RMT for any number of modes N 
in the leads using the Landauer formula, g — 2 Tr(ttf), 
by assuming that the 2N x 2N scattering matrix S = 
r t'\ 

, ) is a random unitary matrix reflecting the er- 
godicity of chaotic scattering M. For single- mode leads, 
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N = 1, this calculation gives P(g) = \{g/2)- 1+ PI 2 §,§, 
shown as dashed lines in Fig. [l[ Notice that the (3 = 1 
distribution is skewed toward smaller conductance, with 
average conductance (g)/3=i = 2/3, while the (3 = 2 dis- 
tribution is constant between and 2 with {g)f3 = 2 = 1. 
The lower average conductance for (3 = 1 results from 
coherent backscattering, analogous to weak localization, 
at B = 0. 

Dephasing, or the loss of quantum coherence, can be 
modeled within RMT by expanding the scattering matrix 
S to include a fictitious voltage lead that supports a num- 
ber of modes 7 V = 2nh/ (t v A), where A = 2irh 2 /m* Adot 
is the spin-degenerate mean level spacing and t v is the 
characteristic dephasing time jl0]|ll]. A recent improve- 
ment to the voltage-probe model that accounts for the 
spatially distributed nature of the dephasing process con- 
siders the limit of a voltage lead supporting an infi- 
nite number of modes, each with vanishing transmis- 
sion, allowing a continuous value for the dimension- 
less dephasing rate 7 V fl2fl . As described below, the 
general effect of dephasing is to make P(g) narrower 
and roughly gaussian, and to reduce the difference in 
mean conductance upon breaking time reversal symme- 
try, 5g = {g)p=2 - (g)f3=i- 




g (e z /h) 



FIG. 1. Conductance distributions for B = (open cir- 
cles), 40 mT (filled circles) and 60 mT (filled triangles) for the 
0.5 fim 2 device at 100 mK, and for B = (open squares) and 
25 mT (filled squares) for the 1.0 /mi 2 device at 45 mK, along 
with the theoretical distributions for fcT/A = 0.61, "f v = 0.3 
(solid curves) and T = 0, 7^ = (dashed curves). Upper left 
inset: Pattern of gates defining each quantum dot. Upper 
right inset: Conductance through 0.5 /1m 2 dot as a function 
of the two shape-distorting gates V g i and V g 2- 

Measurements on two lateral quantum dot with areas 
0.5 /mr 2 (A = 14 meV) and 1.0 /jm 2 (A = 7.1 meV) are 
reported. The devices (see Fig. 1 inset) are defined using 
Cr/Au depletion gates 90 nm above a two-dimensional 
electron gas (2DEG) formed at a GaAs/Alo.3Gao.7As 



heterointerface. Multiple gates are used to allow inde- 
pendent control of the two point-contact leads as well as 
dot shape. Sheet density n = 2 x 10 11 cm -2 and mobil- 
ity fi = 1.4 x 10 5 cm 2 /Vs give an elastic mean free path 
of ~ 1.5 /xm, larger than all device dimensions, so that 
transport within the dots is ballistic. The dots were mea- 
sured in a dilution refrigerator over a range of electron 
temperatures from 45 mK to 750 mK using standard 4- 
wire lock-in techniques at 43 Hz (13 Hz) and less than 7 
/xV (2 /xV) bias voltage for the 0.5 /im 2 (1.0 /im 2 ) device. 
Electron temperature was determined from a fit to the 
average Coulomb blockade peak width, as described in 
& 

Experimental conductance distributions P{g) for 
\B\ « O /A dot {(3 = 1) and \B\ » <f> /A dot {(3 = 2) 
are shown in Fig. [j]. Each histogram contains ~ 1000 
independent samples measured at fixed B drawn from 
the random landscape of conductance fluctuations in the 
space of V g i and V g 2 (inset Fig. [j]). Note that the conduc- 
tance landscape appears random and the average roughly 
constant over the shape-distortion landscape. The asym- 
metry of the (3 = 1 distribution is striking, in contrast 
with previous measurements ]T^JT^ | which found roughly 
gaussian distributions due to thermal averaging and de- 



phasing. The average conductance at (3 = 2 is 
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as expected from RMT, with a small (4%) deviation at 
the lowest temperatures possibly due to imperfect quan- 
tum point contacts enchanced by incipient charging ef- 
fects |l9],^0| . Once time-reveral symmetry is broken, the 
distribution becomes insensitive to magnetic field for the 
relatively small fields used in the experiment (the cy- 
clotron radius is always larger than the dot size). For 
instance, though g(B) are uncorrelated at B = 40 mT 
and 60 mT, P(g) at these magnetic field values are nearly 
identical (<fi /Adot is ^ 8 mT for the 0.5 /xm 2 device). 

Before the experimental distributions can be compared 
to theory it is necessary to measure the dephasing rate, 
since P(g) depends on 7^. Values for ~/ v are measured 
from the weak localization correction to the average con- 
ductivity, dg = {g) [j=2 — (g) p=i using the results of 
Ref. p2[ , as shown in Fig. ||(b). The values for j v (T) 
agree with previous measurements found using a vari- 
ety of magnetotransport methods including weak local- 
ization and the power spectrum of conductance fluctua- 
tions f2l|| , and extend those measurements to lower tem- 
peratures. Unfortunately, while (g) depends on temper- 
ature only through 7^ (which is why weak localization 
is particularly useful for measuring dephasing), the full 
distributions P(g) depend on temperature both implic- 
itly through dephasing and explicitly through thermal 
averaging. The combined effects of dephasing and ther- 
mal smearing must in general be evaluated numerically, 
which we do as follows. A set of values y = J^i Wi(T)xi 
is generated by summing independent samples scj drawn 
from the known distribution P(x,j v ) Jl2fl, weighted by 
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the derivative of the Fermi function, w% = Af'([i + 6] A), 
where /'(e) = 4-(l + e e ' kT )~ 1 , 5 is a binning offset, and 
A is the level broadening, itself dependent on 7 V as de- 
scribed below. By sampling over ensembles of x values, 
a distribution P(y) is obtained (the result is insensitive 
to the choice of 6 for sufficiently large T) . Note that nei- 
ther fluctuations in level spacing nor fluctuations in the 
coupling between the levels and modes in the leads are 
included in this simple model. 

Both dcphasing and temperature averaging tend to 
make P(g) roughly gaussian, in which case P(g) can 
be characterized by its mean and variance. As dis- 
cussed above, the mean, (g), is not affected by ther- 
mal averaging. The variance is reduced both by de- 
phasing, well approximated by the interpolation formula 
Var x = (a + bj v )~ 2 where a = V3 (-y/45/16) and 
b = 1 (y/T/3) for $ = 2 (1) go), and by thermal aver- 
aging, Vary = ^2iWf Vara;. At temperatures exceeding 
the level broadening this sum can be well approximated 
by an integral, 



Vary = A 



Var x = 



6kT 



Var x. 
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The integral form differs from the sum by less than 1% 
for kT > 0.6A. For the devices studied in this paper, this 
condition is easily satisfied, and Eq. (1) is applicable at 
all measured temperatures. 

The thermal averaging procedure given above takes en- 
ergy intervals of size A to be statistically independent. 
Dephasing contributes to level broadening by an amount 
proportional to the dephasing rate, which can be taken 
into account by defining the level broadening to be 



A = A (1+7^/2). 
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Inserting this definition into Eq. (1) reproduces a num- 
ber of previously obtained results for Var g in various 
limits: ( 7v « 1, T « A) @§ ; { lv « 1, T » A) [§; 
(y v > 1, T < A) (1010; and ( 7v > 1, T > A) §. The 
measured variances of the conductance distributions for 
(3 = 1 and (3 — 2 as a function of temperature are com- 
pared to our thermal averaging model in Fig. 2(c). The 
two are in good overall agreement, however, the ratio of 
variances, Var g^i/Var gp=2 shows significant disagree- 
ment between experiment and theory which remains un- 
explained. In particular, the experimental ratio of vari- 
ances is considerably larger than predicted, as seen in 
Fig. §(d). This ratio is an interesting quantity because, 
like 5g, it does not suffer thermal averaging within the 
simple model considered here. Despite the disagreement 
in the ratio of variances, the RMT results for P(g) are 
generally in very good agreement with experiment across 
a broad range of temperatures, as seen in Fig. [|. 

We have also investigated distributions of paramet- 
ric derivatives of conductance with respect to magnetic 
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FIG. 2. (a) Average conductance (g) as a function of tem- 
perature T for B = and 40 mT for the 0.5 fim 2 dot. (b) 
Normalized dephasing rate f v and dephasing time t v deter- 
mined from Sg(T) for the 0.5 fim 2 dot. Not e ag reement with 
previously measured 7 ¥ ,(T) in a 0.4 /j,m 2 dot Ell. (c) Variance 
of conductance for B — and B = 40 mT, corresponding to 
expected variance for f3 — 1 (dashed) and /3 = 2 (solid) in- 
cluding thermal smearing and dephasing effects (see text) . (d) 
Variance ratio Var <j/3=i/Var gp=2 as a function of dephasing 
rate j v . 



field dg/dB and gate voltage, dg/dV g . These quantities 
are of considerable interest as they are the open-system 
analogs of the well-studied "level velocities" dE/dX of 
the energy levels in closed quantum chaotic systems [^J. 
Distributions of parametric derivatives of conductance 
have recently been investigated theoretically using an 
RPA "charged fluid" model, which gave interesting non- 
gaussian distributions for both /? = 1 and (3 — 2 [^2| . 
We have investigated both P(dg/dP>) and P(dg/dV g ) in 
single-mode dots, and find both distributions are well 
described by gaussians due to dephasing and thermal av- 
eraging. Here we focus on P(dg/dV g ), shown in Fig. 4 
for zero and nonzero magnetic field, at T = 0.1 mK and 
0.2 mK. Thermal averaging dominates the width of the 
distribution. We find Vary cx T~ 1,75 . Distributions of 
dg/dB are also roughly gaussian, with variances of 9.65 
mT" 2 at T = 0.1 K and 0.939 mT" 2 at T = 0.5 K. To 
observe deviations from gaussian distributions of para- 
metric derivatives, lower temperatures than those needed 
to see nongaussian P(g) are required. 

In summary, we have measured distribution of conduc- 
tance and its derivatives using large ensembles of shape- 
deformable chaotic GaAs quantum dots with single-mode 
leads. Accurate control of lead conductances and low 
temperature allow the nontrivial predictions of RMT to 
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FIG. 3. Measured conductance distributions for B — 
(open circles) and 40 mT (closed circles) for the 0.5 /im 2 dot 
at four temperatures. Curves show theoretical f3 = 1 (dot- 
ted) and (5 = 2 (dashed) distributions for T = and for 
measurement temperatures (solid). T = distributions are 
determined using the method of Ref. jl2|. Thermal distri- 
butions are calculated according to the sampling procedure 
described in the text. 
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FIG. 4. Conductance derivatives P(dg / dV g ) at B = and 
40 mT for T = 0.1 K and T = 0.2 K. Upper right inset 
shows ~Vax(dg / dV g ) as a function of temperature, which ex- 
hibits a T" 1 ' 75 dependence. Upper left inset shows fluctuating 
conductance g and derivative P(dg/dV g ) as function of gate 
voltage. 



be observed. We find that a thermally averaged RMT 
provides a good description of the measured distribu- 
tions, though some features, in particular the ratio of 
variances Var <7/3=i/Var 3/5=2, are inconsistent with the 
present model and remain to be resolved. 
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